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ABSTRACT 

In the applications like medical radiography enhancing movie features and observing the 
planets it is necessary to enhance the contrast and sharpness of an image. We propose a generalized 
unsharp masking algorithm using the exploratory data model as a unified framework. The proposed 
algorithm is designed to address three issues: 1) simultaneously enhancing contrast and sharpness by 
means of individual treatment of the model component and the residual,2)reducing the halo effect by 
means of an edge-preserving filter, and 3)solving the out of range problem by means of log ratio and 
tangent operations. We present a new system called the tangent system which is based upon a specific 
bregman divergence. Experimental results show that the proposed algorithm is able to significantly 
improve the contrast and sharpness of an image. Using this algorithm user can adjust the two 
parameters the contrast and sharpness to have desired output 

INDEX TERMS: Bregman divergence, exploratory data model, generalized linear system, image 
enhancement, unsharp masking. 



I. INTRODUCTION 

Now a days in digital image processing applications requires an image that contains detailed 
information. So we need to convert a blurry image and undefined image into a sharper image. So in application 
orientation we need to enhance the contrast and sharpness of the image. While enhancing the qualities of the 
image, it is clear that noise components enhancement is undesirable. In this view the enhancement of undershoot 
and overshoot creates the effect called as halo effect. So ideally our algorithm should only enhance the image 
details. To overcome this problem we need to use filters that are not sensitive to noise and do not smooth sharp 
edges. In this paper we use log ratio approach to overcome the out of range problem. In this paper we proposed 
a generalized system which is used for general addition and multiplication which is shown in below figure. 
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And 

a @ x = 0" 1 [a (x)] (2) 

Where x and y are sample signals, a is a real scalar and is non-linear function. And the remarkable property of 
log ratio approach is that it is operable in the region (0, 1) gray scale and can overcome the out of range 
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problem. The following are the main issues in contrast enhancement and sharpness enhancement in current 
existing systems. The existing systems dealt these processes as two tasks which will increase the complexity. 

The contrast enhancement does not lead to the sharpness enhancement. Many of the 
existing systems facing the problem of halo effect. While enhancing the sharpness of an image, in parallel the 
noise of the image also is increased. Though all the systems have enhanced the sharpness and contrast of the 
image, this will give the best result only after careful rescaling process. This was not there in existing system. 
These issues can be solved through our approach using exploratory data model and log ratio approach. 



II. EXPLORATORY DATA ANALYSIS MODEL FOR IMAGE ENHANCEMENT 

A well known idea in exploratory data analysis is to decompose a signal into two parts. From This 
point of view, the output of the filtering process, denoted y = f(x), can be regarded as the part of the image that 
fits the Model. Thus, 

x = y © d (3) 
Where d is called as detail signal (the residual) and defined as d = x Q y where 9 is generalized subtraction. 
The general form of unshaped masking algorithm is as follows. 

V = h(y)©g(d) ~ ~ (4) 
Where v is the output and h(y) and g(d) are the linear or non-linear functions. Explicitly we can say that the 
image model that is being sharpened is the residual. In addition, this model permits the incorporation of contrast 
enhancement by means of a suitable processing function h(y) such as adaptive histogram equalization. As such, 
the generalized algorithm can enhance the overall contrast and sharpness of the image. 

Fig 2. Generalized unsharp masking algorithm block diagram 
The following is the comparison between classical unsharp masking algorithm and generalized unsharp masking 

algorithm. 
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We address the issue of the halo effect by using an edge -preserving filter-the IMF to generate the signal. The 
choice of the IMF is due to its relative simplicity and well studied properties such as the root signals. We 
address the issue of the need for a careful rescaling process by using new operations defined according to the 
log-ratio and new generalized linear system. Since the gray scale set is closed under these new operations We 
address the issue of contrast enhancement and sharpening by using two different processes. The image y is 
processed by adaptive histogram equalization and the output is called h(y). The detail image is processed by 
g(d)= /(oOOtiwhere 7(d) is the adaptive gain and is a function of the amplitude of the detail signal d. The final 
output of the algorithm is then given by 

v = h(y) © [rid} © d] (5) 

III. LOG-RATIO, GENERALIZED LINEAR SYSTEMS AND BREGMAN DIVERGENCE 

The new generalized operations will be defined using the equations (1) and (2). Here these operations 
are defined based on the view vector space used logarithmic image processing. 
A. definitions and properties of log-ratio operations: 

Nonlinear Function: We consider the pixel gray scale xG (0, 1) of an image. For an N-bit image, we can first add 
a very small positive constant to the pixel gray value then scale it by 2" N such that it is in the range (0, 1). The 
nonlinear function is defined as follows: 

0U) = log " (6) 
To simplify notation, we define the ratio of the negative image to the original image as follows: 
X=pG:)=^ (7) 
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Using equation (1) the addition of two gray scales Xj and x 2 can be defined as 
Xl ® Xl = tl^Iw^aOd] = 1^2 (8) 

And the multiplication of the gray scale x by a real scalar or ( — oa < a < +00} is defined using (2) as follows 



a @ x : 



L+jf n 



(9) 



n — 



(ii) 



This operation is called as scalar multiplication and we define a new zero scale denoted by e as follows 
e®x = x ' (10) 

The value of can be defined as follows. 
1/2 <x < 1 
< x < 1/2 
Negative image and subtraction operation: 

The negative of the gray scale, denoted by, is obtained by solving 
x©x=l/2 (12) 

Now we can define the subtraction operation using the addition operation defined in (8) as follows. 
Xl x 2 =x,®(9x 2 ) 



(13) 



Properties: 

Now when positive value a is added to the gray scale x then there will be fluctuations in output image y based 

on the values of a given below. 

Order of reflections of the log-ratio addition. 
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Order of reflections of the log-ratio multiplication. 
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In case of log-ratio addition, based on the value of the constant 'a' the variation of output gray value y is shown 
in the above table. The variations in y depends on the positive and negative nature of the gray scale x and the 
constant 'a'. While coming to log-ratio multiplication the result will be based on the gain value a and the 
polarity of the gray scale x. 
Computation: 

Computations can be directly performed using the new operations. For example, for any real numbers a L 
and H : , the weighted summation is given by 

(«x © *i)®(g a ®*z (14) 
the generalized weighted averaging operation is defined as 

y = (ffij ® *i) © ( «3 ® *3 ) © ( «3 ® 

Where G = (11 1^) 1 * 1 and G =( 11(1 — ^n) c ^) 1/N are the weighted geometric means of the original and the 
negative images, respectively. 

An indirect computational method is through the nonlinear function (6). 

y = _1 [0[Ch 1 ® x L ~) © ( cr, x£ © ( Kfl * n ]} 

(16) 

B. Log-Ratio, the Generalized Linear System and the Bregman Divergence: 
We study the connection between the log-ratio and the Bregman divergence. 
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1) Log-Ratio and the Bregman divergence: 

The Bregman divergence of two vectors x and y, denoted D F (Xll Y), is defined as follows: 

D F (X II Y) = F(X) - F(Y) - (X-Y) T W(Y) (17) 
Where F:K -* Eis a strictly convex and differentiable function defined over an open convex domain K and 
VF is the gradient of F evaluated at point y. the weighted left-sided centroid is given by 
Q. = argmin CE K S„ = i.a n D F (C || X n ) 

= VF- L (^ =1 « n VF(^] (18) 
Comparing equations (17) and (19) we can see that x p is scalar and we can conclude as follows. 
F(x) = J" &tx) dx 

= -xlog(x)-(l-x)log(l-x) (19) 
Where the constant of the indefinite integral is omitted. The previous function F(x) is called the bit entropy and 
the corresponding Bregman divergence is defined as 

D F b ly) - -x log (*) - (1 - x) log(l - x)/(l - >0 (20) 

This is called as logistic loss. Therefore, the log-ratio has an intrinsic connection with the Bregman divergence 
through the generalized weighted average. This connection reveals a geometrical property of the log-ratio which 
uses a particular Bregman divergence to measure the generalized distance between two points. This can be 
compared with the classical weighted average which uses the Euclidean distance. In terms of the loss function, 
while the log-ratio approach uses the logistic loss function, the classical weighted average uses the squared loss 
function. 

2) Generalized Linear Systems and the Bregman Divergence: 

we can derive the connections between the bregman divergence with other established generalized linear 
systems such as the MHS with &(s) = log (x) where x E (0,ocQ and the LIP model &(x~) = -log (1 - x~) 
where x E (— oo, 1). the corresponding bregman divergences are the kullback-Leibler(KL) for MHS 

D t U.y)=xlogi^-U-y) (21) 
And the LIP 

D F ix, y) = (1 - x)log[il - - y>] - [(l - x) - (1 - y}] (22) 

Respectively. 

3) A new generalized system: 

We can define a new generalized system by letting 0GO = VF(i). This approach of defining a generalized 
linear system is based upon using the Bregman divergence as a measure of the distance of two signal samples. 
The measure can be related to the geometrical properties of the signal space. A new generalized linear system 
for solving the out-of -range problem can be developed by using the following Bregman divergence 

D *k,y) = 3r-vT^ ( 23 ) 



Which is generated by the convex function Fix} — —^'(1 — x 2 ~) whose domain is (1,-1). The nonlinear unction 
for E>(jO the corresponding generalized linear system is as follows: 
0U)=™=^= ^ (24) 

In this paper the generalized system is called as tangent system and the scalar addition and multiplication 
(equation (1) and (2)) is called as tangent operations. 

In image processing applications the pixel values from the interval [0,2 N ) is mapped into (1,-1). then the image is 
processed using tangent operations and then the image is mapped to [0,2 N ) as reversible operation. The total 
signal set is confined to (-1,1) and hence using this the out-of -range problem can overcome using log-ratio. We 
can define the negative image and the subtraction operation, and study the order relations for the tangent 
operations. 

IV. PROPOSED ALGORITHM 

A. Dealing with Color Images 

We first convert a color image from the RGB color space to the HSI or the LAB color space. The 
chrominance components such as the H and S components are not processed. After the luminance component is 
processed, the inverse conversion is performed. An enhanced color image in its RGB color space is obtained. 
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The rationale for only processing the luminance component is to avoid a potential problem of altering the white 
balance of the image when the RGB components are processed individually. 

B. Enhancement of the Detail Signal 

1) The Root Signal and the Detail Signal: Let us denote the Median filtering operation as a function y = f(x] 
which maps the input to the output. An IMF operation can be represented as: y^+i = f(yiJ where 
k = 0,1,2,3 is the iteration index andy D = x. The signal y n is usually called the root signal of the filtering 
process if y n +i. = Yn- It is convenient to define the root signal as follows. 

n = min k, subject to H( y^y^-n.) < S 

Where H( yk-Tk+i.} is the suitable measure between the two images and 5 is a user defined threshold. 
It can be easily seen that the definition of the root signal depends upon the threshold. 

2) Adaptive Gain Control: We can see from Fig. 3 that to enhance the detail signal the gain must be greater 
than one. Using a universal gain for the whole image does not lead to good results, because to enhance the small 
details a relatively large gain is required. However, a large gain can lead to the saturation of the detailed signal 
whose values are larger than a certain threshold. Saturation is undesirable because different amplitudes of the 
detail signal are mapped to the same amplitude of either 1 or O.This leads to loss of information. Therefore, the 
gain must be adaptively controlled. In the following, we only describe the gain control algorithm for using with 
the log-ratio operations. Similar algorithm can be easily developed for using with the tangent operations. To 
Control the gain; we first perform a linear mapping of the detail signal d to a new signal c 

c = 2d - 1 (30) 

Such that the dynamic range of c is (-1, 1). A simple idea is to set the gain as a function of the signal c and to 
gradually decrease the gain from its maximum value Ymax when Ic < T to its minimum Yms value when Id— >1. 
More specifically, we propose the following adaptive gain control function 
y(c> = a + /?expt-|c| r (31) 

Where r) is a parameter that controls the rate of decreasing. The two parameters a and P are obtained by solving 
the equations: 

7 (0) = >'\fA?: and y (1) = Ymin ■ For a fixed fj, we can easily determine the two parameters as follows: 
A =fr«« -/]**)/(! (32) 

And 

« = rMAx-2 (33) 
Although both Ymaz an d Ymis could be chosen based upon each individual image processing task, in general it 
is reasonable to set Ymin = 1- This setting follows the intuition that when the amplitude of the detailed signal is 
large enough, it does not need further amplification. For example, we can see that 

^\d\„-+lY® d = lim M^Hi +r ii-ti/<$ = 1 (34) 

As such, the scalar multiplication has little effect. We now study the effect of y and Ymax by setting Ymih = L 

C. Contrast Enhancement of the Root Signal 

For contrast enhancement, we use adaptive histogram equalization implemented by a Matlab function in the 
Image Processing Toolbox. The function, called "adapthisteq," has a parameter controlling the contrast. This 
parameter is determined by the user through experiments to obtain the most visually pleasing result. In our 
simulations, we use default values for other parameters of the function. 

V. RESULTS AND COMPARISON 

We first show the effects of the two contributing parts: contrast enhancement and detail enhancement. 
Contrast enhancement by adaptive Histogram equalization does remove the haze-like effect of the original 
image and contrast of the cloud is also greatly enhanced. Next, we study the impact of the shape of the filter 
mask of the median filter. For comparison we also show the result of replacing the median filter with a linear 
filter having a uniform mask. As we can observe from these results, the use of a linear filter leads to the halo 
effect which appears as a bright line surrounding the relatively dark mountains (for Using a median filter, the 
halo effect is mostly avoided, although for the square and diagonal cross mask there are still a number of spots 
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with very mild halo effects. However, the result from the horizontal-vertical cross mask is almost free of any 
halo effect. In order to completely remove the halo effect, adaptive filter mask selection could be implemented: 
the horizontal- vertical cross mask for strong vertical/horizontal edge, the diagonal cross mask for strong 
diagonal edge and the square mask for the rest of the image. However, in practical application, it may be 
sufficient to use a fixed mask for the whole image to reduce the computational time. We have also performed 
experiments by replacing the log ratio operations with the tangent operations and keeping the same parameter 
settings. We observed that there is no visually significant difference between the results. 



i.iiiiriijal signal -t original signal x 




The above shown graphs are the main differnces between the classical and generalized unsharp masking 
algorithm. Here we show the clear difference between the addition and generalized addition. And at last we 
show that the enhancement is more in generalized unsharp masking algorithm rather than in classical unsharp 
masking algorithm. 
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The above are the two results which processed using generalized unsharp masking algorithm. These two output 
images show that the sharpness and contrast get enhanced in an exponential manner. No rescaling process is 
used here. And the out of range problem is not encountered here since the range of gray scale is (0, 1). 
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